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Abstract 

In this paper, we use a definition of the Debye screening mass in SU{2) 
pure gauge theory at high temperature in 3 + 1 dimensions, that makes use of 
the waUs between different Z2 domains. We show, that this definition enables 
\^ • us to unambiguously obtain the Debye Screening mass from the curvature of 

■ the effective potential for the Polyakov loops calculated in the modified axial 

gauge. Using an approach on asymmetric lattices, we obtain numerical values 
for the Debye mass and compare them to perturbative as well as Monte-Carlo 
results. Our resulting Debye-Masses are significantly higher than the one- 
loop perturbative one. They are within the errors in agreement with those 



obtained by Monte-Carlo methods. Open questions are discussed. 
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1 Introduction 



H , 

5r , The Debye Screening mass in high Temperature SU (2) pure gauge theory has been 

subject of continued interest. The perturbative value in one loop perturbation theory 
is known from calculations of the vacuum polarization tensor as well as other 
observables. 

The perturbative definition of screening phenomena however remains unclear with 
regard to infrared singularities in higher order perturbation theory as well as the 
convergence of the perturbation series in non-abelian gauge theories in general. 
Numerical investigations show, that there is a significant discrepancy in the 
size of the Debye-Mass between calculations on a lattice using Monte-Carlo methods 
and the perturbative result. Up to now, it remains unclear, to what extent higher 
order calculations may overcome these differences or which non-perturbative effects 
might play a role in the process of screening in non-abelian gauge theories even at 
leading order. 

It has been emphasized, that the Haar-Measure associated with the gauge group 
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SU(2) plays a major role when calculating properties of the Polyakov loop in the 
modified axial gauge OqAq = 0. This importance of the Haar measure stems 
from the fact, that in the perturbative limit g ^ the non-abelian nature of the 
Gauge group is inevitably lost, SU{2) [/(l)^^^, thus effectively neglecting effects 
of the non-trivial Haar-measure. Since it is this Haar-Measure, which leads to non- 
trivial topological effects and thus to a feature, which is believed to be specific for a 
non-abehan Gauge theory, it is desirable to take into account its effects also in the 
high-temperature phenomenon of Debye- Screening. 

This paper is organized as follows. After the Introduction in part one, we review a 
definition of the Debye-Mass via the profile of the domain wall between different Z2 
domains of the Polyakov-Loop, which fulfills the requirements of ||^ and suits our 
approach well. Furthermore, we show its consistency with perturbative calculations 
in one loop perturbation theory. In the third part, we summarize some results of 
1^ and § on the method of asymmetric lattices, which will be used further on. 
The fourth part contains our numerical calculations of the curvature of the effective 
Potential for the Polyakov-Loop and the results in terms of the physical temperature. 
The fifth and last part is devoted to a summary of our results and a discussion of 
open problems. 

2 Debye-Mass from the effective Potential of the 
Polyakov-Loops 

There are various suggestions for a definition of the Debye Screening Mass, ranging 
from the exponential falloff of the correlator between two Polyakov loops to the 
static limit of the longitudinal vacuum polarization tensor. A requirement that has 
to be fulfilled by any of them to make sense is its property of being odd under time 
reflection 0. The definition, that we want to use here has been suggested by 
and fulflls the abovementioned requirements. 

It uses the fact, that the center symmetry is broken in the high temperature phase 
in a SU (2) gauge theory and the Polyakov loop may take on one of the two possible 
Z2-values. If there exist two different Z2 phases in space, then a domain wall p{z) 
will develop between them. They can be realized on a lattice by introducing a twist 
in the action and can thus serve as a laboratory for measuring the Debye-Mass. 
The Debye-Mass is then given by the exponential falloff of the profile of the wall 
between different Z2 domains, where here we will investigate only domain walls 
developping in the z direction p{z) 

niD = lim 

The physics of this gauge-invariant non-perturbative definition is clear in the abelian 
case. In the following, we will show in one loop perturbation theory, that it coincides 
with the usual perturbative result in SU{2). 



dzP 
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This goal is achieved by calculating the explicit solution in a Weiss-potential IK 
extending in the z-direction and connecting two different Z2-domains. 
Since in the modified axial gauge the vector potential Uq ~ we have to deal with 
a simple (aj^^-Theory, the field equation of which for a soliton is given by 



dU 
dan 



{z) = 



(2) 



for the potential U[aQ{z)] = ^ — ^^^^^ + const. The solution of this equation is 

aQ{z) = ±mtanh [m{z — zq)] . (3) 



The potential energy of this solution is 

V[ao] = j dx Q(Vao)' + ^(ao) 



(4) 



Using a functional Taylor expansion for small fluctuations r]{z) around the minima 
(ao)o of the potential 



V[ao] = V[iao)o] + I dz- \ 



d^U \ 



viz) + ... 



(5) 



one can see immediately, that the exponential falloff m of the solitonic solution 
is given by the square root of the curvature of the potential at the minimum 



mo 
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Calculating this property for the effective potential in one-loop perturbation theory 



10], one gets the result 



rriD 



(7) 



which coincides with the one-loop result for the Debye-Mass as calculated with other 
definitions. 

It is a peculiarity of the perturbative calculation, that the effect of the nontriv- 
ial Haar-measure is cancelled exactly by the contribution of the longitudinal glu- 
ons, thus leaving the range of integration in the Path integral for the unclear. 
Studies of effective actions allowing for field configurations with a nontrivial space- 
(as compared to the calculation of 
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usmg ao[x) 



const.) 



dependence pl| |ll2| 
have encountered severe difficulties at these points ao = 0, nvr. 
In the following, we will therefore concentrate on the non-perturbative lattice- 
calculation of an effective potential for the Polyakov loops, extracting the Debye- 
mass from its curvature at the minimum subsequently. 
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3 The method of asymmetric lattices 



Although huge efforts have been done to calculate an effective potential for the 
Polyakov loops beyond perturbative methods, no satisfying success has been reached 
yet except for the lattice approach. In this part, we will use the calculation of the 
effective action on an asymmetric lattice as described in [0. 

The essential idea of this method is the fact, that on a highly asymmetric lattice 
with A''^ <C Ns the lattice coupling (3 is replaced by two different spacelike and 

timelike couplings P = p(3s = — , where p is the asymmetry parameter p = \/^. 

p y ps 

The quantum corrections [|13| to these equations will not be considered here. Using 
the abovementioned equations and the fact, that on a highly asymmetric lattice 
/3j /3s, a perturbative series of the action 



x,i \ i i<j / 



in the spacelike coupling jSg can be performed. Only considering the order 0{l3g) for 
the moment, the spacelike link-variables Ui{x,t) of the field strengths originating at 
X in direction i,j respectively 



Goi{x,t) = V{x,t)Ui{x,t + l)V\x + i,t)Ul{x,t), 
Gij{x,t) = Ui{x,t)Uj{x + i,t + l)Uj{x + j,t)UUx,t), 



(9) 



can be integrated. One thus obtains 



/n 


DUx,i;y ^ 


x,i;t 





n p+E 



Nt 



"2 

with Xj the characters in jth representation of the Polyakov loops 

_ , 3^^- 



(10) 



fill 



This result, first obtained by [0] cannot be simplified analytically for general Nt- 
This is why we used 6^ = O^+i to calculate the effective potential 



:i2i 
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for general Nt. [| The 0(/5°) also coincides with the result given in for a strong 
coupling approximation. The reason for this to be the case is the fact that the 0{(3g) 
effectively amounts to neglecting the magnetic parts of the field strength. This is 
exactly how the result 



141 has been obtained. 
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Figure 1: Potential K// in approximation for different couplings b with a 

timelike extension A''^ of the lattice. For increasing coupling b (increasing temper- 
ature), the potential shows a second order phase transition to a phase, where the 
center symmetry is spontaneously broken. 



4 Numerical results 

We now proceed investigating the properties of the effective potential numerically. 
We obtain the effective potential for the Polyakov Loops numerically by only con- 
sidering spatially constant field configurations = 6x+i- 

mean field analysis as performed in for the critical coupling yields the same qualitative 
phase structure but is less practical for calculating the curvature of the potential at the minimum. 
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Then we search for the minima of the effective potential as a function of Pt for dif- 
ferent values of Nt. After that, we calculate the curvature at the minimum by twice 
taking the derivative of Vg// with respect to the argument 9 at the minimum. Next, 
we relate the lattice coupling to the physical temperature ^ and fit the data to the 
functional form of the leading order (one loop) perturbative result for the Debye 
Mass 

mo = C^g{T)T. (13) 
This procedure deserves a more detailed explanation. 

We are eventually interested in the prefactor C. Previous lattice calculations 0, 0] 
suggested a C G [1.6, 2] as compared to the one loop perturbative result C = 1. 
For being able to make reliable predictions about the size of this prefactor, one has 
to compare results at the same physical temperature. We follow here an approach 
outlined in [|T^. 

Violations of asymptotic scaling are taken into account by considering the general 
Ansatz 

aA = R{g')\{g') (14) 



487r2' 3 V167r2/ 

The function \{gt) parametrizes the asymptotic scaling violations, the values for 6o 
and bi are found perturbatively. For X{gt), we use an exponential ansatz 

\{g^) = exp (:^id^g^ + d^gf + d^g!)] . (17) 



With T = Tf—, we obtain 



2bl 



1 ^^2^T, 



KglH- (18) 



The value of the coupling g^^ at the critical Temperature of the deconfinement 
phase transition can be obtained numerically by inspection of the effective potentials 
for various A''^ (see figure |I]). We thus arrive at the fit parameters di = 0.0068, d2 = 
-0.001173, 4 = 6.4818 10"^ in equation 

Now we can extract the temperature in units of the critical temperature for a given 
Nt. It is given by 

T _ R{gl)X{gl) 
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Nt 




mo 


^TTS 


1 
i 


n QQ 
u,oy 


1 QS 

i,yo 


i,Uio 


2 


1,29 


1,81 


0,889 


3 


2,09 


1,80 


0,983 


4 


2,54 


1,64 


1,156 


5 


3,18 


1,64 


1,330 


8 


5,37 


1,48 


1,556 



Table 1: Results 



This sets us in the position to relate the Temperature T to the coupling [3t. This 
temperature scale can now be used to relate the lattice Debye mass given by the 
curvature of the effective potential ([T2|) at the minimum in terms of the asymmetric 
lattice coupling f3t to the physical temperature T. 

We now use the functional dependence of the Debye-mass on g{T) and T as known 
from one- loop perturbation theory to determine C in equation (|l^). To this end, we 
need the coupling as a function of the temperature g{T). For this running coupling, 
we use the two-loop formula 



9{T) 



11 
12^ 



In ( 27r-^ 
T. A 



17 
44^ 



In 



21n( 27r— ^ 
A 



(20) 



Using (|Tj), (|19D and (|20D in (p!3|), we can determine the prefactor C which gives the 
quotient of m£, with the perturbative value m^*'. 

Carrying out this procedure for several values of Nt in the temperature range 2 < 
^ < 9>,we arrive at the values for the prefactor C in the one loop perturbative Debye 
mass given in table |l|. These results are shown in figure El. The values for that 
we obtained with our fit are also indicated in the table. A comparison with the 

A 



numerical result = 19.8231 [|r^ has been performed. Thus the deviation from 
one gives the degree of accuracy and shows, that it agrees well with the expection 
for small lattices, whereas it strongly disagrees for larger lattices. This will also be 
adressed in the next section. 



5 Discussion and Outlook 

Using a definition of the Debye-Mass which allows to extract it from the gauge in- 
variant effective potential for the Polyakov loops, we showed by explicit solution, 
that the perturbative result is reproduced, when the Weiss-Potential is considered. 
We then used an asymmetric lattice to determine the effective potential nonper- 
turbatively. To this end, we only considered the leading order contribution in an 
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Figure 2: Comparison of the Debye masses obtained for various values of Nt. The 
Value of and the older one of as well as the perturbative value are indicated, 
too. 



expansion in the space-like lattice coupling. This approximation is believed to be of 
good accuracy for highly asymmetric lattices. This can also be shown numerically 
by inspection of the curvatures of the corresponding effective potentials in 0{(3^) 
and 0{l3l) for different values of the asymmetry parameter p (which ultimately give 
the Debye mass). The result is depicted in figure 0. The influence of the next to 
leading order on the temperature scale is not taken into account here. 

Our investigation then proceeded by numerical calculation of the Debye-Mass. 
The results of our calculation yield a prefactor which is significantly larger than the 
perturbative result. It is within the error bars due to the fits, that we performed 
compatible with the result suggested by full Monte-Carlo calculations 0, which is 
also shown in figure |^. 

The major shortcoming of the method used here is the fact, that no Monte-Carlo 
values of ^ are available for highly asymmetric lattices (only these can be hoped to 
give reliable predictions for the Debye-Mass). Thus there is no control as to whether 
the temperature scale is of good accuracy and whether the results for C are reliable. 
In addition to that, since we used the 0(/3°) of the asymmetric expansion only, 
we are effectively working in the strong coupling regime and neglecting magnetic 
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Figure 3: Difference of tfie curvature of tlie effective potentials in V^'^l , and 



0{(3g), V^j:j , at the minimum for Nt = 5 (full line) and Nt = 3 (dashed line), 
|: = 11, 5. For p = ^ = 1, the deviation is 22% at A^j = 5 and 24% at Nt = 3, for 
p = 20 only 0,81% and 3,8% respectively. 



terms in the action. Due to this, our results are not obtained in the scaling regime, 
which can also be seen from the funtion A((?^) (|1^, which fits the asymptotic scal- 
ing violation and is not constant. This behaviour of the function parametrizing 
the asymptotic scaling violation is also seen in full numerical investigations [00. 
However, this problem can easily be overcome by incorporating the contribution of 
the 0{(3g) terms. Our investigation of the errors of the truncation of the series in (3s 
shows, that the effect on the Debye Mass will be small on the scale of the discrep- 
ancy between perturbative continuum result and full lattice Monte-Carlo result. 
Furthermore, we had to do several fits, which are of varying accuracy. In this par- 
ticular case, we performed the analysis for rather small values of f3t, because the 
expansion of the action in (3s is better for small (3. In addition to that, we used per- 
turbative results for R{g'^) as well as for g{T). It is unclear, whether the coupling 
range that we investigated is in the perturbative regime already. 
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On the other hand, the resuhs show, that there is a strong influence of the wall at 
PL^ ±1 in the effective Potential for the PL associated to the non-trivial Haar- 
measure on the size of the inverse screening length. Since the fact, that this wall 
does not vanish is a major property of non-perturbative treatments of the action 
(as has also been seen in a path-integral approach in 1 -|- 1 dimensions |^ and in 
a new, remarkable calculation I^H), the difference in size of the Debye-Mass can 



be attributed to this Haar-measure, which inevitably modifies the curvature at the 
minimum of the effective potential. 
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